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^ \ We study the resonances of the quantum kicked rotor subjected to an excitation 

;> that follows an aperiodic Fibonacci prescription. In such a case the secondary reso- 

nances show a sub-ballistic behavior like the quantum walk with the same aperiodic 
prescription for the coin. The principal resonances maintain the well-known ballistic 



t~^. ' behavior. 
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The quantum kicked rotor (QKR) is considered as the paradigm of period- 
ically driven systems in the study of chaos at the quantum level pp. Two 
important characteristics of the behavior of the quantum kicked rotor are dy- 
namical localization (DL) and the ballistic spreading of the resonances [2]. 
These behaviors are quite different and depend on whether the period of the 
kick t is a rational or irrational multiple of 4-7T (in convenient units). For ra- 
tional multiples, the behavior of the system is resonant and has no classical 
analog. For irrational multiples the average energy of the system grows in a 
diffusive manner for a short time and then the diffusion stops and localiza- 
tion appears. From the theoretical point of view the two types of values of r 
determine the spectral properties of the Hamiltonian, for irrational multiples 
the energy spectrum is purely discrete and for rational multiples it contains a 
continuous part. Both quantum resonance and DL can be seen as interference 
phenomena, the first is a constructive interference effect and the second is a 
destructive one. 
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We recently developed [3] a new generalized discrete time quantum walk (QW) 
on the line and showed that this model has a dynamical behavior analogous 
to that of a QKR: depending on the values of a parameter there are either 
quantum resonances or DL. This modified QW has been mapped into a one- 
dimensional Anderson model [I], as was previously done in the case of the 
QKR [5]. For some exceptional values of the parameter, which correspond to 
resonant behavior, the standard deviation grows linearly in time. We showed 
that the usual discrete time QW on the line becomes a particular case of 
resonance of the modified QW. 

The concept of QW on the line is a subject that has drawn much attention in 
quantum computing [6J. It has been introduced in 1993 [7|I8] . as counterpart 
of the classical random walk. A classical random walk is defined in terms 
of the probabilities for a particle to do a step to the left or to the right 
but the QW is described in terms of probability amplitudes. Many classical 
algorithms are based on classical random walks, then it is possible that future 
quantum algorithms will be based on the quantum random walk. One of the 
most striking properties of the QW is its ability to spread over the line linearly 
in time as characterized by the standard deviation cr(t) ~ t, while its classical 
analog spreads out as the square root of time (cr(t) ~ t 1 / 2 ). 

In this work we establish another aspect of the parallelism between the QKR 
and the QW following the proposal of Ribeiro et al. |9j. In [9] the quantum 
coin operator of the QW is arranged in aperiodic sequences using the Fi- 
bonacci prescription and this leads to a sub-ballistic wave function spreading 
(cr(t) ~ t c , 1/2 < c < 1). We show that the secondary resonances of the QKR, 
excited with the same Fibonacci prescription, have also sub-ballistic behav- 
ior like the usual QW, but the primary resonances maintain the well-known 
ballistic behavior. Then the parallelism with the usual QW on the line is re- 
stricted to the secondary resonances of the QKR. On the other hand Casati et 
al. [10] studied the dynamics of the QKR also kicked according to a Fibonacci 
sequence, but outside the resonant regime, and they found sub-diffusive be- 
havior for small kicking strengths; in this sense our work can be considered as 
complementary to theirs. 

The QKR Hamiltonian is 

p2 oo 

H = — + Kcos6j25(t-nT) (1) 

2-* n =l 

where the external kicks occur at times t = nT with n integer and T the 
kick period, I is the moment of inertia of the rotor, P the angular momen- 
tum operator, K the strength parameter and 9 the angular position. In the 
angular momentum representation, P\£) = £h\£), the wave-vector is = 
IX-oo^00l^) and the average energy is E(t) = {9f\H\9) = e EZ-oo ^ M*)| 2 , 
where e = h 2 /2I. Using the Schrodinger equation the quantum map is readily 
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obtained from the Hamiltonian ([I]) 



(2) 



j=-oo 



where the matrix element of the time step evolution operator U (k) is 

U ej = i-^e- ij2£T / H J j - t (K) 1 (3) 

J m is the mth order cylindrical Bessel function and the argument is the di- 
mensionless kick strength k = Kj h. The resonance condition does not depend 
on k and takes place when the frequency of the driving force is commensu- 
rable with the frequencies of the free rotor. Inspection of eq.© shows that the 
resonant values of the scale parameter r = eT/2h are the set of the rational 
multiples of Air, t = 4tt p/q. In what follows we assume, that the resonance 
condition is satisfied, therefore the evolution operator depends on k, p and q. 
We call a resonance primary when p/q is an integer and secondary when it is 
not. 
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Fig. 1. Standard deviation a(t) as a function of time, a) random sequence for the pri- 
mary resonance p/q = 1; b) periodical sequence for the primary resonance p/q = 1; 
c) Fibonacci sequence for the secondary resonance p/q = 1/3; d) random sequence 
for the secondary resonance p/q = 1/3. 

With the aim to generate the dynamics of the system we consider two values 
of the strength parameter k, K\ and k 2 and combine the corresponding time 
step operators U («i) and U (k 2 ) in a large sequence that takes three different 
forms, namely periodic, random and quasi-periodic. In this way we generate 
three types of unitary evolution operators. With these operators we compute, 
for several thousands of T, the wave function spreading as measured by the 

exponent c in a(t) = JEZ-oc? 2 M*)? ~ t C - 
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Fig. 2. The exponent c as a function of the ratio p/q that identifies the secondary 
resonance. The probability distribution was calculated, for the Fibonacci quasi-pe- 
riodic sequence, with k\ = 5, K2 = 10. The value p/q = 1/2 does not appear because 
it corresponds to an antiresonance. The symmetry of the figure is explained in the 
text. 

We make a numerical study of the long time behavior of the parameter c for the 
three types of sequences. In the periodic case, for different sequences, the bal- 
listic behavior, c ~ 1, is still found for all resonance conditions. In the random 
case we performed several calculation where the number of operators U («i) 
and U (/C2) in the random sequence were taken in different ratios, obtaining in 
always c ~ 1 for the primary resonances and c ~ 1/2 for the secondary ones. 
The quasi-periodic case was performed using the Fibonacci prescription given 
in [9] where U n = \J n -\U n -2 with Uq = U (ki) and U\ = U (1C2); for example, 
applying the above rule six times we get UiU UiUiU UiU Ui and this sequence 
gives the dynamical evolution up to t = 8T. In this case, two types of results 
have been obtained: ballistic behavior (c ~ 1) in the primary resonances and 
sub-ballistic behavior (1/2 < c < 1) in the secondary resonances. The stan- 
dard deviation is plotted in Fig. [H which displays the qualitative differences 
between the periodic, random and quasi-periodic cases. 

It is possible to prove analytically why the ballistic behavior is maintained 
in the primary resonances for all three types of sequences. For the primary 
resonance with the initial condition a\ (0) = 5iq the solution of the map eq.([2]), 
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with a given sequence of U («) operators, is az(t n ) = (—i) J? n(j) I = 

\j=o J 

(— i) Je(rrii <%i + m 2 k 2 ), where takes the values K\ or k 2 involved in 
the sequence and m\ and m 2 are the numbers of occurrences of K\ and k 2 
respectively, mi+ m 2 = n. With the above amplitudes, the standard deviation 
is obtained easily as 

VyiT) = n (4) 

where a = m\jn and /3 = m 2 /n are the relative weights of kicks «i and k 2 
respectively. Thus c = 1 and, additionally, we have obtained explicitly the 
ballistic diffusion coefficient D = a{nT)/n = (azKi + (3n 2 ) /v2. 

The QKR has an exceptional behavior when p/q = 1/2, called antiresonance, 
its characteristic being that the system returns to the initial state every two 
periods. In this case, given a sequence of step operators, the solution of the 

map eq.©, is ae(t n ) = {—i) 1 Jt \ {"^Y K {j) ) an d the standard deviation 



is a{nT) = {.—i) K (j)/v2- Therefore, for all sequences, the antiresonance 

in the case n\ ^ k 2 has a similar behavior as the primary resonances. It is im- 
portant to emphasize, that the conceptual key for the ballistic behavior of the 
primary resonances (and the antiresonance too) resides in the commutativity 
between operators U (ki) and U(k 2 ). This is not the case for the secondary 
resonances where the commutator does not vanish. 



The sub-ballistic behavior takes place for all secondary resonances in the quasi- 
periodic case. In Fig. [2] the exponent c is plotted as a function of the ratio 
p/q that identifies the secondary resonance. From this figure we can conclude 
that the exponent c depends on both p and q, but there is a trivial symmetry 
in the time step evolution operator eq.([3]) when p/q is changed by (q — p)/q, 
this is the reason why e.g. the value of c for p/q = 1/5 is the same as for 
p/q = 4/5. We have also found that the dependence of c on the values of Ki 
and k 2 is not smooth in general. The condition i%i = k 2 corresponds to the 
periodic case and the expected ballistic behavior is obtained. Fig. [3] shows the 
cut of the surface c(ki,k 2 ) with the plane K\ = —k 2 for p/q = 1/3, making 
clearly apparent the sub-ballistic behavior of the system for all k. We have 
also studied higher moments of order four and six. The asymptotic behavior 
of these moments is consistent with the power-law behavior of the second 
moment, i.e. all the moments obtained with the Fibonacci prescription have 
smaller exponents than those obtained with a periodical sequence. 

At this point it is possible to ask oneself if this sub-ballistic behavior has 
a quantum origin or if it also appears in the classical world. The classical 
equations of motion for the Hamiltonian eq. (CQ) are given by the standard map 
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Fig. 3. The exponent c, that characterizes the power law of the standard deviation, 
as a function of the strength parameter k of the kicked rotator. The dynamical 
evolution is obtained by the operators U(ki) and U(k2) following a Fibonacci suc- 
cession. 



where n indicates the time step. In the phase space (9 n ,P n ) of the standard map 
a beautiful diagram is obtained where the Kolmogorov-Arnold-Moser (KAM) 
tori prevent diffusion in angular momentum for small values of K < K cri . 
The classical model equivalent to the quantum model developed in this work 
is obtained when the strength parameter K takes the values K\ and K 2 in a 
given sequence in correlation with the time step. We have worked numerically 
with the map eq.(jSJ) in the periodic, random and aperiodic (following the 
Fibonacci prescription) sequences. In the periodic case the KAM surfaces are 
established, but in the random and aperiodic cases the KAM surfaces are 
broken and the classical diffusion, is established. Then we can conclude that 
the sub-ballistic behavior of the model is a quantum phenomenon. 

Quantum resonances and DL, that were at first established in numerical and 
theoretical form, have been experimentally observed more than ten years ago 
in samples of cold atoms interacting with a far-detuned standing wave of laser 



P n+ i = P n + K sin 6 n 
9 n+ i = 6 n + P n+ i mod 27r 



(5) 
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light [T2lll3|ll4j . This notable series of experiments have drawn much attention 
because they may be establishing both the conceptual and experimental basis 
of quantum computers [Toll 1 611 17|TT5] . Experimentally only the primary reso- 
nances of the QKR are easily observable, but recently, Kanem et al. [IS], have 
observed secondary resonances. On the other hand, several experiments have 
been proposed [20|21||22||23] to construct models of QWs. These proposals, at 
first sight, have some common elements with the experimental implementation 
of the QKR in the recent past. Thus, in this theoretical and experimental frame 
where the QKR and the QW have equivalent behaviors and their experimen- 
tal facilities have many elements in common, the question to be posed is: may 
the QKR be considered as the fundamental model of QW? More experimental 
and theoretical work is still necessary in order to answer this question. 

In summary, we have developed an aperiodic QKR model and established 
a deeper equivalence between the QKR and QW. In this model we found a 
new sub-ballistic behavior in the secondary resonances, a(t) ~ t c with 1/2 < 
c < 1, where c depends on k, p, q. The usual ballistic behavior of the QKR 
model is retained in the primary resonances, then these resonances are robust 
because they are not affected by the type of sequences of the operators U (^i) 
and U (^2), periodic, random or aperiodic. We explained this robustness by 
the commutativity of the U operators and we obtained the ballistic diffusion 
coefficient analytically. Finally, we showed numerically that the sub-ballistic 
behavior is a quantum interference phenomenon that has no classical analogue. 

We acknowledge the support from PEDECIBA and PDT S/C/IF/54/5 and 
V.M. acknowledges the support of D.F.P.D.-ANEP. 
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